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ABSTRACT 

We present a new geometrical method for measuring the distance to the 
Galactic center (Rq) by solving for the Keplerian orbit of individual stars bound 
to the black hole associated with the Sgr A* from radial velocity and proper 
motion measurements. We identify three stars to which the method may be 
applied, and show that 1-5% accuracy of R can be expected after 15 years of 
observing, and 0.5-2% after 30 years of observing, depending on what the orbital 
parameters of these three stars turn out to be. Combining the measurements 
of the three stars with favorable orbital parameters leads to even more precise 
values. In the example that we present, such combined solution yields 4% 
accuracy already by the year 2002. All these estimates assume that annual 
position measurements will continue to be made with the a p = 2 mas precision 
recently reported by Ghez et al. The precision of the distance measurement is 
relatively insensitive to the radial velocity errors, provided that the latter are 
less than 50kms~ 1 . Besides potentially giving an estimate of Ro that is better 
than any currently in use, the greatest advantage of this method is that it is free 
from systematic errors. 

Subject headings: celestial mechanics, stellar dynamics — binaries: 
visual — Galaxy: center — Galaxy: fundamental parameters 
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1. Introduction 

The distance to the center of the Galaxy (Rq) is to Galactic astronomy what the Hubble 
constant (H ) is to extragalacatic astronomy and cosmology. Many of the parameters we 
determine for the Galaxy depend on R . Fixing its value will directly give us a precise value 
of circular velocity 60 from the measured proper motion of Sgr A*. The circular velocity 
and distance in turn enter in estimates of the mass of the Galaxy. Having a precise value of 
Ro would allow us to recalibrate distance indicators such as RR Lyraes, and compare them 
with local RR Lyraes. This would give us an idea of applicability of the local RR Lyrae 
calibration to globular clusters and the LMC. This can to some extent also be done for the 
Cepheids, which are an important step in the extragalactic distance ladder. 

As with the Hubble constant, the distance to the center of the Galaxy is still debated 
(although with perhaps less controversy), 70 years after it was first estimated by Shapley 
(1918). In 1993, Reid presented an overview of then current methods and suggested a 
weighted average of Rq = 8.0 ± 0.5 kpc. This is a convenient value which is still used by 
many researchers (although the IAU maintains its recommended value at Rq = 8.5 kpc). In 
the closing sentence of his review, Reid (1993) expects Rq to be known with 3% uncertainty 
by the end of the decade. How close are we to this goal today? In the period after the 
Reid's (1993) article there were four new estimates of R , one of which was based on a new 
method. 

The first new estimate uses RR Lyrae observations in the near-infrared (IR). Apparent 
magnitudes suffer less interstellar extinction in the near-IR than in the optical. Smaller 
brightness amplitudes in the near-IR also permit an estimate of the mean magnitude with 
fewer measurements. Using 58 RR Lyraes in Baade's Window, Carney et al. (1995) report 
Ro = 7.8 ± 0.4 kpc which is consistent with Reid's value and has 5% uncertainty. However, 
the real uncertainties turn out to be greater, as this result depends on the distance-scale 
model used. In fact, the authors report two possible results that differ by 14%. Metzger, 
Caldwell & Schechter (1998) measured Rq based on Cepheids and a model of the Galactic 
rotation. They find Rq = 7.66 ± 0.32 kpc. The 4% uncertainty is a combination of many 
different systematic errors (Cepheid calibration, reddening, etc) that enter in the final 
result. The value of Rq itself is also model dependent, and an additional uncertainty comes 
from an unknown component of the ellipticity of the stellar orbits. 

The Cepheid distance scale was recalibrated using the Hipparcos satellite. This led 
to a new estimate of Rq by Feast & Whitelock (1998). They report a "long" scale of 
Rq = 8.5 ± 0.5 kpc, i.e., a 6% uncertainty. 

Hipparcos data also facilitated the introduction of a new method using red clump stars. 
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Paczyriski & Stanek (1998) used Hipparcos parallaxes to calibrate nearby red clump stars 
and applied this calibration to find the distances of the red clump stars near the Galactic 
center. Due to the great number of these stars in Baade's Window, the formal uncertainty 
of Rq is very small. However, various systematics produce a final estimate of Rq = 8.4 ± 0.4 
kpc, i.e. 5% accuracy. Later, this result was corrected to Ro = 8.2 ±0.15 ±0.15 kpc (Stanek 
& Garnavich 1998). 

As we can see, there is still no agreement on a single value of Rq, and all of the 
estimates have significant uncertainties. Notably, none of these methods falls into a category 
of primary methods - those that are not using a "standard candle" or a Galactic rotation 
model (Reid 1993). One would prefer to use a direct method in order to avoid systematics. 
Even the H 2 maser method (Genzel et al. 1981; Reid 1993) is model dependent. Also, the 
accuracy of this method seems to be limited to about 15%. 

An attractive way to directly determine Ro is to measure the trigonometric parallax of 
Sgr A* using VLBI. However, significant technical obstacles remain for this approach (Reid 
1996, private communication). The Space Interferometry Mission (SIM) will be sensitive 
enough to measure parallaxes with 3% precision at the Galactocentric distance, but since 
it operates in the optical part of the spectrum, SIM cannot look directly at the Galactic 
center. 

Here we suggest a new, purely geometric method, which became feasible with the 
recent advances in speckle imaging of the stars surrounding the Galactic center (Eckart & 
Genzel 1997; Genzel et al. 1997; Ghez et al. 1998). In essence, it is the classic visual binary 
method which has been used for decades to measure the masses of and distances to binary 
stars. We apply this method to stars orbiting around the massive object in the center of 
the Galaxy (presumably a black hole). This method requires observations of the proper 
motions and radial velocities of several stars near Sgr A* during the next decade or so, in 
order to constrain the shape and the physical size of the orbit. These, together with the 
angular size of the orbit will in turn yield R . This is not a statistical method: each star 
individually leads to a value of R whose uncertainties depend only on the precision of the 
proper motions and radial velocity measurements. Its basic advantages are that it is free of 
systematic errors and that it can possibly measure Ro accurate to 1-5% after 15 years of 
observing, and 0.5-2% after 30 years of observing. 

The orbital solution to a visual binary automatically gives the heliocentric radial 
velocity of the center of mass which in this case is nearly identical to Sgr A*. Since the 
velocity of the Local Standard of Rest (LSR) is known, the radial velocity of Sgr A* 
immediately gives the velocity of the LSR relative to Sgr A*. If this velocity is found to 
differ significantly from zero, it can only mean that either Sgr A* is not at rest with respect 



-4- 



to the Galaxy or that the LSR is moving radially with respect to the Galaxy. Gould & 
Popowski (1998) find that the LSR is moving radially outward relative to local halo stars 
at —4.0 ± 8.5 kms" 1 . Measurement of the velocity of Sgr A* would therefore allow one to 
distinguish between these two possibilities. 



2. Physical Principles 

Eckart & Genzel (1997) and Ghez et al. (1998) have measured the proper motions of 
many stars within 1" of Sgr A*. These proper motions are actually Keplerian orbital motion 
of stars around the central mass of M = 2.6 x 10 6 M Q (this assumption will be discussed 
in more detail in § 5.). As yet, the time baseline of these observations (~ 2 years) is not 
sufficient to show these proper motions as making an arc around Sgr A*. Let us consider 
star S0-1 from the Ghez et al. (1998) catalogue (SI in Genzel et al. 1997 catalog), at a 
separation of 0."1 from Sgr A*, the closest so far detected. Of all the stars in the vicinity 
of Sgr A*, S0-1 is the most likely to show the signs of orbital motion in near-term future 
observations. How long do we have to wait until this happens? If we assume a circular 
orbit, and that the projected separation at which we now see S0-1 is close to the angular 
radius of its orbit, then Kepler's law gives us a period of 17 years. This is just a rough 
estimate, but it gives us an idea that the full orbit of S0-1 can be traced in a reasonable 
amount of time. 

What will learning about the orbit of S0-1 tell us about RqI For simplicity, assume that 
S0-1 is in a circular orbit of period P = 15 yr, with a major-axis a = RqOl, where a is the 
current separation of a = 100 mas. Let the angle of inclination be i, so that the projected 
semi-minor axis is b = a cost. The observed dispersion in radial velocity measurements is 
related to these quantities by 

1/2 2 2 1/2 27rasku 27rasin? 
[var^d)] 1 / 2 = ((vl d ) - (v Tad ) 2 )^ 2 = — = R (1) 

Let us for the moment assume that P and var(f ra( j) can be measured with perfect accuracy. 
Then the fractional error in the determination of R® is equal to the fractional error in the 
measurements of asini. The fractional error in a is simply 

5a <j„ /2 yr „ / a a \ ( a \ _1 ' '" 



0.007 



a a V P \2mas/ VlOOmas/ \15y r 



implying that the method is potentially very precise. However, for low-inclination orbits, 
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the fractional error in ctsini is much greater than in a. Specifically 

S(asmi) 5(Ja' 2 — (3 2 ) J a 2 + 8 2 . VT+ cos 2 i 5 a , oX 

— ^— = — / 2 m = ~2 — w 6a = — ' 3 

a sm i y/a 2 — p 2 cr — p z sm % a 

where (3 = b/R . That is, the fractional error in R is larger than the fractional error in a 
by a factor of 2.5 for % = 45°, and by a factor of 5.3 for % = 30°. The fractional error in the 
period will be 5P/P ~ 5a/ a, so it is appropriate to ignore it in this simplified analysis. 

How good do the radial velocity measurements have to be so that they do not dominate 
the error in the determination of The fractional error in the dispersion measurement 
must be smaller than the fractional error in asini. That is 

5 v /var(w rad ) _ ( j^ ad (P/2yr)- 1 /2 ^(q; sin _ y/1 + cos 2 i 5a 



^/var(v rad ) v osini «sini sin i 



a 



where v = 2na/P, and where we have assumed that there are annual radial velocity 
measurements with accuracy <J Viad - Equation (4) can be expressed 



o p VT+ cos 2 i _j (o p /a\ V esc 2 i + cot 2 i ( v \ , , 

o v , < — v = 56kms — - — — — . (5) 

a sini V0.02; 2 V1400kms~V KJ 

Hence, the velocity errors will not dominate unless they are greater than 50kms _1 . 

This simple analysis for the case of a circular orbit demonstrates the potency of the 
method. In the next section we will do the exact (although intuitively less clear) treatment 
of this problem in the case of orbits with arbitrary eccentricity and orientation. We present 
the accuracies obtainable for R from the three stars closest to Sgr A*. In § 5., we discuss 
the validity of our basic assumptions and some specific details. 



3. The Method 

The distance to Sgr A*, R , is one of the 13 parameters required to describe the 
Star-Sgr A* binary orbit, and the precision of its measurements depends not only on the 
observational errors, but on the actual values of the 13 parameters, i.e. on the specific 
orbit. We can only specify the precision of this distance measurement as a function of these 
13 parameters. At first sight, this appears unwieldy. However, as we now show, 11 of the 
parameters are known sufficiently well that the uncertainty in R does not depend on them 
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in a major way. They can therefore be considered fixed, and it is only necessary to specify 
the uncertainty in R as a function of the two remaining free parameters. 

The 13 parameters can be taken to be: (1) the six phase-space co-ordinates of Sgr A*, 
(2) the six phase-space co-ordinates of the orbiting star relative to Sgr A*, and (3) the 
mass M of Sgr A*. The six phase-space co-ordinates of Sgr A* and the mass of Sgr A* 
can be considered known. For example, we take the distance to be 8 kpc. However, if the 
distance is 5% higher, 8.4 kpc, our estimate of the distance error will be off by a similar 
amount, which is quite small compared to the factor ~ 5 range of values we will explore for 
the distance error. Similarly, we take the mass of Sgr A* to be M = 2.6 x 1O 6 M , but if 
it turns out to be 5% higher or lower, our estimate of the distance error will have to be 
rescaled by a similar amount. Nothing in the error analysis depends in any way on the 
values of the two-dimensional angular position, the two-dimensional angular velocity and 
the radial velocity of Sgr A*, so these quantities can be truly ignored. 

Four quantities are measured: the two-dimensional angular position and velocity of the 
star. Unfortunately, these measurements are relative to the frame of the near-IR speckle 
images, not Sgr A*. However, Menten et al. (1997) have aligned this IR frame with the radio 
frame in which Sgr A*'s position and proper motion are measured. Using this alignment, 
Ghez et al. (1998) fix the position of Sgr A* to ±10 mas. Ultimately, the visual-binary 
orbital solution will fix this position with an order of magnitude better precision. 

For the orbit to be completely specified, two additional parameters are required. For 
convenience, we take these two parameters to be the period P and inclination i. Thus, for 
each P and % we construct an orbit that is consistent with the observable quantities: the 
two-dimensional angular position and two-dimensional angular velocity. Not all values in 
the P — i parameter space produce an orbit. When they do there are four possible solutions. 
There are two sets of two solutions that each comprise a pair of orbits that are symmetrical 
in the plane of the sky and so have the same error structure. Thus, there are effectively two 
different solutions (two different orbits) for each P and i. We calculate orbital parameters 
for all allowed orbits having periods P < 40 yr (we choose this as an upper limit in our 
study since longer periods will take too long to constrain the orbit observationally) . 

After we find the parameters for all the possible orbits of a given star we simulate 
a set of observations made once a year. Observations consist of measuring the position 
and the radial velocity. We introduce the radial velocity measurements in the fifth year 
of observations (i.e. in the year 2000, since we set t obs = at the epoch of the first 
Ghez et al. (1998) observations, i.e., the year 1995), when we expect adaptive-optics 
spectroscopy to become good enough to produce these data. We assume the current 
positional measurement accuracy (2 mas) throughout, and an error in the radial velocity 
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measurements of <7„ rad = 20kms 1 . The choice of a Viad will be discussed in the next section. 

We determine the evolution of the accuracy of -R as the observations accumulate by 
calculating the 13 x 13 covariance matrix for the n = 13 parameters of the model. To do so, 
we first write the three observables F m (t) (to = 1, 2, 3 for F m = p X) p y , t> rad ) as functions of 
the time t and of the 13 parameters which we label a i: F m (t; ai, . . . , ai 3 ). We then linearize 
F m in the neighborhood of the solutions, a*. 

n 

Funit) = £(* - OM) + F m (t; a*,,..., a* 3 ) (6) 
i=i 

where 

Qpm 

We then specify the solution using seven classical parameters for the orbit of the star 
relative to Sgr A*, and six to specify the three-dimensional position and velocity of Sgr A* 
itself. Specifically, we have 

F m (t) = G m (t) + p m , + fi m t, (to = 1, 2) (8) 

F\t) = R Q d -^ + v rfl (9) 

where (pi,P2)o is the initial position of Sgr A* and (/xi,/x 2 ) is its proper motion relative 
to the IR star frame, t> r ,o is the heliocentric radial velocity of Sgr A*, and G m (t) is the 
three-dimensional orbit of the star relative to Sgr A*, in angular units. That is, 



G m (t) = £^(*i)^(fe)rcj.(*3)ir(f), (io) 

pqr 

H\t) =a(cos^-e), H 2 (t) = a{\ - e 2 ) 1/2 sin^, H 3 (t) = 0, (11) 

ip — e sin ip = uot + 0. (12) 



Here, a is the angular semi-major axis, e is the eccentricity, u = 2n/P, <p is the phase of 
the orbit at t — 0, and 1Z 1 , 1Z 2 , and 1Z 3 are the rotations which we take to be around the 3, 
1, and 3 axes respectively. 
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We then create the matrices 6, and c with elements 

Qpm Qpm ^ 
da i 9a 3 {°kY 

where a™ is the error in the k-th measurement of F m (tk)- The index k steps through the 
"data" - in our case the observations made each year. The square of fractional uncertainty 
of the parameter ctj is just the co variance matrix element Cjj. One of these parameters is 
Rq. This allows us to predict the accuracy of finding R after t obs years of observing. In the 
next section we present the results of these calculations for the three stars closest to Sgr A*. 



4. Predictions of the Uncertainty of R 

How does the accuracy of the R determination evolve as the observations accumulate 
at the rate of one per year? Let us illustrate this using the star SO-1, presently the closest 
to Sgr A*. For each P and % we will for now consider only one of the two possible orbits. 

In Figure 1 we fix the inclination at % — 45° and explore all the possible orbits with this 
inclination. Different lines correspond to different possible periods in integer increments. 
The short-dashed line corresponds to the shortest possible period of P = 20 yr, while the 
long-dashed one represents the orbit of P = 40 yr, the upper limit that we explore. Each 
line represents the evolution of the fractional uncertainty of R , as the observing time t Q b s 
progresses. The first thing to notice is that for the orbits of a fixed inclination, <jr /Ro 
will not depend strongly on the period. The plot shows that 5% accuracy is achieved after 
about 12 years of observing, 2% after some 20 years, and 1% accuracy after approximately 
30 years. After this time, which corresponds to observing a star throughout its entire orbit, 

— 1/2 

a plateau is reached, i.e. additional observations reduce the uncertainty of R only as t obs . 

We now explore the set of possible orbits with % = 55°. This is close to the highest 
possible inclination of i max = 57° that star SO-1 could have. The evolution of (Jrq/Rq is 
shown in Figure 2. The shortest possible orbit (short-dashed line) now has P = 32 yr. 
Compared to i — 45° case, the uncertainties of R now depend even more weakly on period. 
The trend of a Ri) /R evolution is also different. An accuracy of 5% is reached after 16 
years, i.e. later than in the case of % — 45, while 2% and 1% accuracies are obtained after 
about the same time as for % = 45°, 20 and 30 years, respectively. While in the case of 
% = 45° orbits, the maximum accuracy is reached at 0.7% level, now it is much lower, 0.4%. 
This illustrates that this method can potentially yield a very low uncertainty for R : it all 
depends on what the orbit of this star (and the other two) turn out to be. It should be 
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noted that the second set of solutions (for a given P and i) essentially give the same trends 
of a Ro /R evolution. 

As we mentioned above, the radial- velocity errors are assumed to be a Vvad = 20kms~ 1 . 
Since these measurements have not yet been made, we do not know what their precision 
is going to be. However, we can explore how sensitive the determination of R® is to a 
range of different values of o"t, rad . In Figure 3 we show the evolution of ct Ro /Rq for SO-1 
for the P = 30 yr, % = 45° orbit (one of the lines in Figure 1) for various assumed values 
cr Viad = 5, 10,20,50 and 100 kms -1 . The basic result is that cr Ro /Ro is not strongly affected 
by <7„ rad over this range. At late times (t b s > 30 yr) <j Ro /Rq depends only very weakly on 
<J Vrad . In the region of t Q b s where this dependence is strongest we see that a Vigd = 20kms _1 
(thick line) is the highest value that is still close to the lines of low a Viad . That is why we 
adopt a Vigd = 20kms _1 throughout this paper. The <J Ro /Ro dependence on cr v ^ d is much 
less for the orbits with short periods (regardless of inclination), and only somewhat greater 
for longer periods. 

Finally, we move on to exploring the entire P — i parameter space. We do so by making 
contour plots for S0-1 in the P — i plane of a Ro /R after some fixed observing time. (We 
use contours showing 10, 5, 2, 1, and 0.5% error). We examine two cases: t Q b s = 15 yr 
and t bs = 30 yr. The outermost contour in all of the contour plots actually limits the 
parameter space of possible orbits. Figure 4 displays contours of ctr /-R f° r orbits with 
different P and i after t Q b s = 15 yr. It shows that orbits of higher inclination allow more 
precise determination of Rq. Since for each P and i we have two possible orbits, in Figure 
5 we show the contours for the second set of solutions, also after t Q b s = 15 yr. The region 
where <J Ro /R < 5% is smaller in this case. After t obs = 30 yr, the two solutions (Figs. 6,7) 
are much more alike. The contours of high R precision (better than 1%) are again in the 
high inclination range and towards shorter periods. This corresponds to short orbits (full 
orbit is traced faster) that are away from face-on orientation where the fractional error in 
asini is high. See § 2. 

Star S0-1 is the closest star from Sgr A* (0."114) found to date and has the highest 
measured two-dimensional velocity (~ 1400 km s -1 ) (Ghez et al. 1998). Therefore, it 
represents the first logical choice for the analysis we showed. However, we find two other 
stars in the Ghez et el. (1998) catalog that can also be used for finding Rq in t Q b s = 30 
yr, and one of them produces useful results in t Q b s = 15 yr. These are SO-2 (0."151 from 
Sgr A*) and SO-3 (0."218 from Sgr A*). They are also known as S2 and S4 in the Genzel 
et al. 1997 catalog. (S3 from Genzel et al. 1997 is another candidate, but its detection is 
ambiguous in the Ghez et al. 1998 catalog which we use throughout this paper. 

We find (but do not show) that the P — i parameter space for the star SO-2 looks 
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substantially different from S0-1. Highly inclined orbits (almost edge-on) are allowed for 
fairly short orbital periods. After t obs = 15 yr some orbits already produce R estimates 
with better that 1% errors. This is several times better than in the case of S0-1. After 
^obs = 30 yr, some orbits of SO-2 even give a better than 0.5% estimate of Rq. Finally, for 
the star SO-3, after t Q b s = 30 yr some orbits yield a better than 1% estimate of Rq. No 
useful estimate of Rq can be obtained from SO-3 after t obs = 15 yr. This is not surprising 
since the shortest possible orbits have P = 19 yr. 

How much would the three stars together, compared to each of them individually, 
contribute to a more certain value of RqI One can simply take a weighted average of the 
estimates of Rq from each single star. To illustrate this, in Figure 8 (note the extent of 
y-axis compared to Figs. 1-3) we choose some values for the orbital parameters for each of 
the stars S0-1 (P = 30 yr, i = 45°), SO-2 (P = 20 yr, i = 70°) and SO-3 (P = 30 yr, i = 75°). 
The evolution of (Jrq/Rq from each of them is shown by thin solid lines. The uncertainty 
obtained by averaging is shown by a thick dashed line. However, observing the three stars, 
even for a short time, can produce an even better estimate of Rq than the simple weighted 
average. This is because partial orbits of more than one star can together constrain the 
position, velocity and mass of Sgr A* much better than each orbit individual: each star 
must orbit a common mass at a common center. We show the reduced uncertainty by the 
solid thick line. In the first few years this combined solution produces a very significant 
improvement compared to the naive statistical treatment. Afterwards, it remains better by 
some 40%. Note that in this example, 4% precision can be obtained already in the eighth 
year of observation, and 1% precision only two years after that, i.e. by 2003 and 2005, 
respectively. Ultimately, 0.2% accuracy of Rq could be achieved. We also performed the 
same analysis excluding the observations of SO-3. We find that SO-3 contributes significantly 
(mostly by constraing the mass of Sgr A*) even at times when the observations of SO-3 
alone do not tell us much about Rq. This indicates that perhaps even the stars farther away, 
that alone cannot determine Rq (those with P > 40 yr), might contribute by constraining 
the mass of Sgr A*. Therefore, it would be advisable to keep track of other stars near Sgr 
A* in addition to the three closest. 

We note that in this example, the radial velocity of Sgr A* is measured to 5 km s _1 
precision after 20 years and 3kms _1 after 30 years, roughly half the measurement errors of 
Gould & Popowski (1998). See § 1. 
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5. Discussion and Conclusions 

Besides the possibility of giving us a precise measurement of Rq, this method is very 
powerful because it does not rely on some model or a calibration. It has a single underlying 
assumption: stars in the vicinity of Sgr A* move in a Keplerian potential. We believe 
that this is a robust assumption which is consistent with the available observational data. 
Nevertheless, it is possible in principle that in addition to the point mass at the center of 
the potential there is a continuous mass distribution composed of many faint stars. This 
would cause the orbits to precess. However, this motion, if significant, would show itself 
after the orbit has been fully traced, and can then easily be accounted for. We would still 
obtain a precise estimate of R . 

Recently, Munyaneza, Tsiklauri & Viollier (1998) have explored an exotic version 
of this scenario with an extended object at the Galactic center composed of degenerate, 
self-gravitating heavy neutrino matter. They show that such an object would alter the 
orbit of SO-1 so as to make it distinctly different from a Keplerian orbit, and that positional 
measurements alone will be able to resolve this over the next 10 years. Distinguishing 
between the black-hole and extended-object scenario will become even easier once the radial 
velocity curve of SO-1 is measured. 

In § 2 we argued that in the case of circular orbits, only relatively crude radial velocity 
measurements (o"i, rad < 50kms _1 ) are required. We now investigate the role of radial 
velocity measurements more closely. In Figure 3 we show the fractional error in Rq as a 
function of time for SO-1 assuming i = 45° and P = 30 yr, and for various values of <J Vrad . 
Also shown in this plot (dashed curve) is the fractional error in the quantity a sin if P. This 
would be exactly equal to ct Ro /Rq, if there were perfect knowledge of var(v rad ) (i.e. perfect 
velocity measurements) but none of the radial-velocity information were used to disentangle 
degeneracies among the orbital prameters. 

At late times (t obs > P) all the curves converge with the dashed curve lying slightly 
below the others. This shows that high-precision velocity measurements are indeed of small 
value if the observations are carried out for at least one period. However, at early times, 
the dashed curve lies well above the others. This shows that for observations shorter than 
the period, the radial velocity measurements play a significant role in disentangling the 
degeneracies among the orbital parameters, and hence can significantly help determine the 
distance Rq. 

The success of the method we present, in terms of precision, will eventually depend on 
how lucky we are to see the stars in the vicinity of Sgr A* in favorable orbits. In any event, 
this method will produce an estimate of the distance to the Galactic center that is free from 
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the various systematic effects that affect other methods currently in use. 
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Fig. 1. — Evolution of Rq fractional uncertainty from SO-1 orbits with % = 45° (t obs is 
in years). Different lines correspond to different integer orbital periods. The short-dashed 
line corresponds to the minimum possible period of P = 20 yr, while the long-dashed line 
represents the orbit of P = 40 yr, our upper limit. 



Fig. 2. — Evolution of R fractional uncertainty from SO-1 orbits with % = 55° (t obs is 
in years). Different lines correspond to different integer orbital periods. The short-dashed 
line corresponds to the minimum possible period of P = 32 yr, while the long-dashed line 
represents the orbit of P = 40 yr, our upper limit. 



Fig. 3. — Evolution of R fractional uncertainty from SO-1 orbit (P = 30 yr, % = 45°, t obs 
is in years) with different assumptions for the measured a Viad = 5, 10,20,50 and 100 km s^ 1 
(upper solid line corresponds to 100 km s -1 ). Bold line represents <7„ rad = 20km s -1 . Dashed 
line is the fractional error in asini /P. See § 5. 



Fig. 4. — Contour plot of R fractional uncertainty from S0-1 orbit, after t obs = 15 yr. 
Solution 1. 



Fig. 5. — Contour plot of Rq fractional uncertainty from S0-1 orbit, after t obs = 15 yr. 
Solution 2. 



Fig. 6. — Contour plot of R fractional uncertainty from S0-1 orbit, after t obs = 30 yr. 
Solution I. 



Fig. 7. — Contour plot of Rq fractional uncertainty from S0-1 orbit, after t obs = 30 yr. 
Solution 2. 



Fig. 8. — Evolution of Rq fractional uncertainty from each of the stars (Pi = 30 yr, i\ = 45°; 
P2 = 20 yr, %i = 70°; P3 = 30 yr, i 3 = 75°) separately (thin solid lines), and when combined 
in a straightforward statistical way (thick dashed line), and a combined solution with a 
common position, velocity and mass of Sgr A* (thick solid line). 
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